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Abstract. The phenomenological models of convection use characteristic length scales they do not determine but that are 
chosen to fit solar or stellar observations. We investigate if changes of these length scales are required between the Sun and low 
mass stars on the red giant branch (RGB). The question is addressed jointly in the frameworks of the mixing length theory and 
of the full spectrum of turbulence model. For both models, the convective length scale is assumed to be a fixed fraction of the 
local pressure scale height. 

We use constraints coming from the observed effective temperatures and linear radii independently. We rely on a sample of 38 
nearby giants and subgiants for which surface temperatures and luminosities are known accurately and the radii are determined 
through interferometry to better than 10%. For the few cases where the stellar masses were determined by asteroseismological 
measurements, we computed dedicated models. First we calibrate the solar models. Then, with the same physics, we compute 
RGB models for masses between 0.9M o and 2.5M Q and metallicities ranging from [Fe/H] = -0.34 to solar. The evolution is 
followed up to 10 3 L . A special attention is given to the opacities and to the non grey atmosphere models used as boundary 
conditions for which the model of convection is the same as in the interior. 

We find that for both the mixing length theory and the full spectrum of turbulence model the characteristic solar length scale for 
convection has to be slightly reduced to fit the lower edge of the observed RGB. The corresponding models also better match the 
expected mass distribution on the RGB and are in better agreement to the seismic constraints. These results are robust whether 
effective temperatures determined spectroscopically or radii determined interferometrically are used. 

Key words. Physical data and processes : convection. Stars individual : low mass - evolution - red giants. Techniques: high 
angular resolution 



1. Introduction 

The modelling of solar and stellar surface convection requires 
hydrodynamical radiative transfer computations (Nordlund et 
al.|2009). These tridimensional calculations are extremely time 
consuming because convection zones motions are character- 
ized by short time scales and enormous turbulence. Therefore, 
in spite of the current computational power, the direct mod- 
elling of convection does not rely on first principles only but 
still necessitates simplifying assumptions. 

Direct observations can constrain the properties of outer 
convection in stars. The Sun's radius is routinely used to cali- 
brate the main free parameter of the mixing length theory, a m i t , 
the characteristic length scale of convection in this model be- 
ing generaly a-mitHp where H p denotes the local pressure scale 
heigth. It is customary to consider that the solar a m \ t applies all 
the way through stellar evolution and some recent studies have 
suggested this to be a valid hypothesis for low-mass red giant 
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stars (Ferraro et al. I2006I VandenBerg et al. [2008). However 
hydrodynamical convection calculations (Ludwig et al. 1999) 
show that slightly different or m i t should be used depending on 
the effective temperature (hereafter T e jj) and the surface grav- 
ity (hereafter log g). Furthermore, the solar a m i t appears in- 
appropriate to describe outer convection in red giants in the 
3 - 20 M mass range (Stothers & Chin ll997i l and during the 
final stages of stellar evolution (D'Antona & Mazzitelli 11 9961 1. 
This is one reason why it is interesting to check what the nearby 
Galactic red giants can tell us about the mixing length the- 
ory (hereafter MLT). Besides the MLT there is another local 
treatment of convection that is increasingly used: the full spec- 
trum of turbulence model (hereafter FST, Canuto, Goldman & 
Mazzitelli 1 19961 1 . This approach is more physically consistent 
than the MLT and shows better agreement to numerous sets of 
observations (Mazzitelli 1 1 9991 and references therein, Samadi 
etal.|2006]l. 

As red giants and solar surface conditions (T e ff, log g) 
strongly differ, they offer an adequate opportunity to check 
the MLT and the FST empirically. Such analyses were made 
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in the recent years. However in the case of the nearby red gi- 
ants, the usual observational constraints, absolute luminosity 
and T e ff are now supplemented by direct interferometric radii 
measurements and the new asterosesimic constraints. In this 
work we aim at using these new data to estimate a m \ t and a cgm 
and their possible variation from the Sun to Galactic red gi- 
ants. The analysis is made for both convection models using 
constraints on effective temperatures and radii independently. 
In section 15.31 we give the physical inputs to our stellar evo- 
lution code. We describe in detail the treatment of convection 
and the atmosphere modelling. In section[3]we present and dis- 
cuss the observation sample. The calibration of the convective 
length scales relying on the Sun is made in section|4] Section[5] 
examines the MLT and the FST model of convection under the 
constraints coming from the red giant branch stars. We discuss 
our results and conclude in section [6] 

2. The evolution code and the models 

We use the CESAM code (Morel 1997 ), an hydrostatic one di- 
mensional stellar evolution codqj. However the standard ver- 
sion of CESAM is significantly modified here: to allow the 
computation of the stellar structure in the late stages of evo- 
lution we changed the usual integration variables of CESAM 
(see appendix 1). The evolution is computed from the zero age 
main sequence up to a luminosity of 1O 3 L on the red giant 
branch (hereafter RGB) for stars less massive than 2M . For 
more massive stars, the models are evolved to the end of helium 
core burning on the early asymptotic giant branch (hereafter 
eAGB). The nuclear reaction rates are based on the NACRE 
compilation (Angulo et al. 1999). The network includes the 
proton-proton chains, the CNO cycle, and 3a and l2 C(a, y) 16 
reactions. The neutrino losses are computed according to the 
analytical fits of Itoh et al. (119961 ). Microscopic diffusion is 
taken into account in the solar calibration models (® and 
more generally for all the models down to an effective temper- 
ature of 5000 K. This temperature roughly corresponds to the 
maximum extent in mass of the outer convection zone near the 
base of the RGB, the so-called first dredge up. Then because of 
the deep convective mixing, the subsequent evolution of diffu- 
sive and non diffusive models become undistinguishabk^l (see 
Salaris et al. (2002 ) and references herein) and diffusion is not 
taken into account below T e ff = 5000K. 

The calculation of solar radius or RGB/eAGB radii models 
critically depends on the physical assumptions and inputs used 
in the outer layers. There are four of them. 

i) The opacities: The low temperature opacities domi- 
nated by metallicity effects are critical in determining the 
effective temperature of giant stars at a given luminos- 
ity (Salaris et al. 120021 ). Below log T = 3.75 we use 
the data from Ferguson et al. (2005 ) (publicly available at 

1 A list of the scientific publications using CESAM is available at 
http://www.oca.eu/morel/articles.html 

2 However we recall that even if the diffusive and non diffusive 
models reach the same RGB they do not do so at the same age because 
the main sequence downward diffusion of helium and heavy elements 
speeds up the evolution on that stage. 



http://webs.wichita.edu/physics/opacity/). In the higher tem- 
perature regime we use the OPAL opacities computed from 
the Lawrence Livermore National Laboratory web interface 
(http://physci.llnl.gov/Research/OPAL/new.html). Both opac- 
ity sets correspond to the solar metal repartition advocated by 
Asplund et al. (2005). The electron conductive opacities that 
are only significant in the degenerate helium cores of red gi- 
ants follow the Itoh et al. ([1983') prescription. We model stars 
well below the RGB tip and using the newer prescription of 
conductive opacities by Cassisi et al. (2007) would likely leave 
the evolution tracks unchanged as noticed by these authors. 

ii) The convection: Two simplified local treatments of con- 
vection are investigated: the mixing length theory and the full 
spectrum of turbulence model. We chose a prescription of the 
MLT very similar to that of Bohm-Vitense ( 119581 ) and whose 
exact description is given in the appendix of Piau et al. (120051 ) . 
The formulation of the FST is that of Canuto, Goldman & 
Mazzitelli ( 1996 ) and referred as CGM hereafter. The detailed 
equations are provided in the second appendix to this article. 
For the CGM approach, we adopt a characteristic length scale 
A = CcgmHp in the equations -where H p is the pressure scale 
height-. The CGM version we use is a simplified version rather 
than the true version of the CGM where A is the distance to 
the convection/radiation regime boundary. Finally we have two 
length scale parameters for the two treatments of convection: 
ffmit and a cgm . In the occurence of convective cores, i.e. for stars 
more massive than 1 .2M Q , we assume convective overshooting 
ranging from 0.1H p to 0.2H p (Claret |2007]). 

Hi) The atmosphere modelling: We define the atmosphere 
as the region where the Rosseland optical depth is less than 
Tb = 20. The outer boundary conditions to the internal structure 
are taken at this depth where the diffusion equation becomes 
valid (Morel et al. 1 19941 ). In late type stars, the convection 
straddles such an atmosphere/interior boundary and because 
the MLT and the CGM model predict different temperature gra- 
dients, it is important to use the same treatment of convection 
above and below this limit (Montalban et al. 120011 Montalban 
et al. 2004). Besides this, molecular lines have a significant im- 
pact on atmosphere structures when T e ff < 5000 K. For these 
reasons we introduce two series of non grey atmosphere models 
as outer boundary conditions. The first series of temperature- 
optical depth relations (T(r) 4 = Tj ff f gr jd(T)) is computed with 
the PHOENIX/1 D atmosphere code where the convection is 
handled using the MLT. The second series of temperature- 
optical depth relations is computed with the ATLAS 12 atmo- 
sphere code (Castelli 120051 We modified ATLAS 12 in order to 
use the CGM version of the FST model of convection instead 
of the MLT. The change is directly inspired from the imple- 
mentation of the FST to the Atlas9 code by Kupka Q996). The 
characteristic length scales adopted for the convection in both 
atmospheric sets is 0.5H p as suggested by the solar and the cool 
dwarfs Balmer lines (Samadi et al. 2006, Gardiner et al. |1999l ). 
Using the same convection model in the atmosphere and inte- 
rior is necessary for consistency. In the proper treatment of the 
CGM, the convection length scale A is the distance to the con- 
vection/radiation regime boundary. This mimics the increase 
of the convective efficiency with depth (Heiter et al. 2002) and 
the smooth change of the thermal gradient associated. However 



L. Piau et al.: Surface convection and red giants radii measurements 



3 



this treatment is not handled by our code where A = ffcgmHp. 
In this framework, it is impossible to fit the solar radius when 
using in the interior the (small) atmosphere a cgm required by 
the B aimer lines. The sudden change in a cgm (or ormit) from the 
atmosphere to the interior should induce a discontinuity in the 
thermal gradient at the limit between the interior and the atmo- 
sphere. In order to smooth it out we perform a linear interpola- 
tion with optical depth on the temperature gradients: provided 
V a and V; are the temperature gradients (V = computed 
using respectively the atmosphere and interior formalisms, we 
compute the gradient V at the the optical depth r < t\, follow- 
ing: 

V = (l-x(T))V a + x(T)Vi 

where x(t) = max(0, (t — l)/(Tb - 1)). Thus the atmospheric 
gradient is taken where t < 1 , the interior gradient where r > 
Tb and and linear combination of them where 1 < t < tj,. 

Both atmosphere grids were computed with the solar 
Asplund et al. (2005 ) composition. In both of them there is a 
model every 100 K on the [4000,6400] K effective temperature 
range and every 0.5 on the [0.5,5.5] range of decimal logarithm 
of surface gravity (log g). For every stellar model we perform 
a linear interpolation in T e ff and log g using the four closest 
neighbours on the grid in order to get the f(r) theoretical rela- 
tions. Together with T e ff this relation provides the temperature 
atmosphere profile: T(t) 4 = T 4 ff f(r). For the most massive stars 
modelled, T e ff > 6400 K on the main sequence and the f gr ;d(T) 
function is not available. In this case we use the f gr id(r) at 6400 
K and for the closest log g. 

iv) The equation of state (hereafter EOS): We use the OPAL 
2005 EOS tables for Population I stars. Although the effect of 
the EOS is somewhat less important than the effects of opacities 
or atmosphere boundary conditions (Salaris et al. 120021 1 it still 
determines the temperature gradient in the convection zone and 
therefore affects the radius. 

Let us make a final comment on the composition of the 
models before moving to the next section. The atmosphere T(t) 
relations are especially computed for the Asplund et al. (2005) 
composition and metal repartition. The EOS and opacities both 
in low and high temperature regimes assume the same metal 
repartition. However we do not only consider the solar compo- 
sition. The spectroscopic analyses of the red giants of our sam- 
ple suggest them to have on average a slightly lower metallicity 
than the Sun and small but non negligible composition differ- 
ences ( 35.21 ). The total metal mass fraction Z and helium mass 
fraction Y are changed to explore the impact of composition. 

3. The observations 

The selection of our sample of stars was done in two steps. We 
first queried the CHARM2 catalogue (Richichi et al. 120051) to 
obtain all direct measurements of giant and subgiant angular 
diameters up to 2004, with effective temperatures in the range 
5000 to 5500 K. The giants from Hutter et al.'s (fT989l Mark III 
survey were removed from the sample, as most of these mea- 
surements present strong biases due to calibration uncertainties 
(see their Sect. IV for details). We then searched the literature 
for more recent observations, and added the measurement of 



y Sge by Wittkowski et al. (f2006), 5 Eri and £ Hya by Thevenin 
et al. J2005T ) and the recent high accuracy CHARA/FLUOR 
measurements of eOph (Mazumdar et al. 120091) and ^Ser 
(Merand et al. 2009). When several independent angular diam- 
eter measurements were available for the same star, we com- 
bined them into a single value taking into account their orig- 
inal error bars and the consistency of the different measure- 
ments. The conversion of uniform disk angular diameters to 
limb darkened values was done using linear limb darkening co- 
efficients by Claret ( 1995), that are based on stellar atmosphere 
models by Kurucz (Kurucz 1993 ). We are aware that the plane 
parallel ATLAS9 limb darkening coefficients by Claret (1995 ) 
are not optimal for cool, giant stars with low effective gravi- 
ties, but the difference with PHOENIX models is very small for 
the selected sample. This difference amounts typically to a few 
per mille of the stellar size. Choosing the Claret d 19951 ) values 
has the advantage of preserving the internal consistency of the 
sample, as the ATLAS9 models are the standard for the quoted 
interferometric measurements. A comparison of the measured 
angular diameters with predictions from the Fk{ V - K) surface 
brightness-color relations calibrated by Kervella et al. (|2004a) 
on Cepheid supergiants is presented in Fig. [T] The agreement 
is satisfactory within the uncertainties, with no systematic bias 
at a level of 2%. It is interesting to remark that these relations, 
established for supergiants, are mostly identical to the relations 
calibrated by Kervella et al. (2004b ) using dwarf and subgiants. 
This is an indication that the infrared surface brightness-color 
relations are universal for all classes of F-K stars. 

The spectroscopic effective temperature T e ff , effective grav- 
ity logg and metallicity [Fe/H] of each star of our sample 
were taken from Cayrel de Strobel et al.'s ( 120011) catalogue, 
except for ySge and ^Hya for which we used the T e (f val- 
ues determined in the articles reporting their interferomet- 
ric angular diameters. These estimates were originally ob- 
tained by McWilliam ( [19901 ), Thevenin & Idiart ( [1999] ) and 
Mallik (1998). Similarly the T e ff, metallicity and luminosity 
of eOph are adapted from De Ridder et al. (120061) (also in 
agreement with Mazumdar et al. (2009)). The bolometric lu- 
minosity was estimated using the K band photometry and the 
corresponding bolometric corrections BC{K) by Houdashelt et 
al. (120001 ). For ySge, we extrapolated the BC(K) value for 
T eff = 3805 K (Wittkowski et al. 120061 For eOph and 5 Eri, 
we adopted the T e ff metallicity and luminosity from De Ridder 
et al. (2006 ) and Thevenin et al. (2005) respectively. To derive 
the linear radius, we used the parallaxes from the Hipparcos 
catalogue ( I1997I ). We kept only the stars for which the relative 
uncertainty on the linear radius is lower than 10%. The dis- 
tances to the selected stars range from 11 to 1 lOpc. Thanks to 
this proximity, we neglected the interstellar reddening for the 
computation of the bolometric luminosity. 

Our final sample contains 38 giant and subgiant stars with 
spectral types G5 to M0. It is interesting to remark that three 
stars of our sample have known asteroseismic oscillation fre- 
quencies: 5 Eri, £Hya, eOph. This allowed an accurate deter- 
mination of their masses through a combination of their ra- 
dius and the large frequency spacing (Thevenin et al. [2005 ; 
Mazumdar et al. |2009l ), that are estimated respectively to 1.215, 
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Fig. 1. Comparison of the limb darkened angular diameters of 
the stars of our sample from direct interferometric measure- 
ments (horizontal axis) and photometric estimates based on the 
Fk(V - K) infrared surface brightness-color relation (IRSB) 
from Kervella et al. (I2004al >. The lower plot shows the residu- 
als as a function of the measured angular diameter. 



2.65 (no error bars mentioned by the authors) and 1.85 + 
0.05 Mg. 



Table 1. Observational and physical parameters of our selected sample of giant and subgiant stars with interferometrically measured angular diameters. 



Star 


HD 




Sp. Type 


LD 




71 


Radius 




l°gg 


[Fe/H] 


m v 


m K 


bc(a:) 


M bol 














(mas) 




(mas) 


(Re) 




(K) 














(L e ) 




p Her 


HD 


161797 


G5IV 


1.94 ± 


0.03 


119.05 ±0.62 


1.75 i 


:0.03 


5520 


3.70 


0.04 


3.41 


1.74 


1.56 


3.68 ± 0.05 


2.7 ±0.1 




6 En 


HD 


23249 


K0IV 


2.39 ± 


0.03 


110.58 ±0.88 


2.32 ± 


:0.03 


5055 


3.79 


0.13 


3.51 


1.40 


1.84 


3.46 ± 0.05 


3.3 ±0.2 




/JAql 


HD 


188512 


G8IV 


2.23 ± 


0.10 


72.95 ± 0.83 


3.28 ± 


: 0.15 


5100 


3.60 


-0.30 


3.71 


1.65 


1.84 


2.81 ±0.06 


6.0 ± 0.3 




7/Cep 


HD 


198149 


K0IV 


2.67 ± 


0.04 


69.73 ± 0.49 


4.12 ± 


:0.07 


4950 


3.41 


-0.32 


3.41 


1.22 


1.84 


2.28 ± 0.05 


9.7 ± 0.5 




7/Ser 


HD 


168723 


K0III-IV 


2.94 ± 


0.01 


53.93 ±0.18 


5.87 ± 


:0.03 


4890 


3.21 


-0.42 


3.26 


1.05 


1.84 


1.55 ± 0.05 


19 ± 1 




46 Lmi 


HD 


94264 


K0III 


2.56 ± 


0.03 


33.40 ±0.78 


8.22 ± 


:0.22 


4670 


2.96 


-0.20 


3.83 


1.30 


2.00 


0.91 ±0.07 


34 ±2 




cr Gem 


HD 


62044 


K1III 


2.33 ± 


0.05 


26.68 ± 0.79 


9.38 ± 


:0.35 


4500 


2.40 


-0.02 


4.28 


1.74 


2.16 


1.03 ± 0.08 


31 ±2 




eAql 


HD 


176411 


K1III 


2.00 ± 


0.12 


21.22 ±0.77 


10.14 


±0.70 


4760 


2.91 


0.00 


4.03 


1.79 


2.00 


0.42 ± 0.09 


54 ±5 


P 


37Tau 


HD 


25604 


K0III 


1.70 ± 


0.09 


18.04 ± 0.84 


10.15 


±0.69 


4700 


2.77 


0.01 


4.35 


2.03 


2.00 


0.31 ±0.11 


60 ±6 


£Hya 


HD 


100407 


G7III 


2.39 ± 


0.02 


25.23 ± 0.83 


10.18 


±0.35 


5037 


2.83 


-0.04 


3.54 


1.45 


1.84 


0.30 ± 0.09 


60 ±5 


w 


eOph 


HD 


146791 


G9.5IIIb 


2.96 + 


0.01 


30.62 ± 0.20 


10.39 


±0.07 


4955 


2.89 


-0.27 


3.24 


0.98 


1.84 


0.25 ± 0.05 


63 ±3 


a 


(5Ari 


HD 


19787 


K2III 


1.88 ± 


0.13 


19.44 ± 1.23 


10.42 


±0.97 


4810 


2.93 


-0.03 


4.35 


2.17 


2.00 


0.61 ±0.14 


45 ±6 




0CrB 


HD 


136512 


K0III 


1.18 ± 


0.07 


11.90 ±0.74 


10.64 


±0.93 


4730 


2.75 


-0.44 


5.52 


2.93 


2.00 


0.30 ±0.14 


60 ±8 


in 


eCyg 


HD 


197989 


K0III 


4.56 ± 


0.02 


45.26 ±0.53 


10.82 


±0.14 


4730 


2.89 


-0.27 


2.50 


-0.01 


2.00 


0.26 ± 0.06 


62 ±3 


I 


<5Tau 


HD 


27697 


K0III 


2.20 ± 


0.07 


21.29 ±0.93 


11.12 


±0.61 


4810 


2.93 


-0.03 


3.76 


1.64 


2.00 


0.28 ±0.11 


62 ±6 


re 

re 


39ari 


HD 


17361 


K1.5III 


1.89 ± 


0.11 


18.06 ± 0.84 


11.26 


±0.82 


4600 


2.85 


-0.02 


4.51 


2.10 


2.16 


0.54 ±0.11 


48 ±5 


re 
o 


<5'Tau 


HD 


27697 


K0III 


2.34 ± 


0.02 


21.29 ±0.93 


11.81 


±0.53 


4030 


1.83 


-0.29 


3.76 


1.64 


2.51 


0.79 ±0.11 


39 ±4 


3 


12Aql 


HD 


176678 


K1III 


2.42 ± 


0.01 


21.95 ±0.92 


11.85 


±0.50 


4600 


2.75 


-0.19 


4.03 


1.47 


2.16 


0.34 ±0.10 


58 ±6 


re 


6 And 


HD 


3627 


K3III 


4.09 ± 


0.02 


32.19 ±0.68 


13.64 


±0.30 


4360 


2.37 


0.04 


3.28 


0.47 


2.16 


0.17 ±0.07 


68 ±4 


! 


a Ari 


HD 


12929 


K2III 


6.84 ± 


0.02 


49.48 ± 0.99 


14.85 


±0.30 


4480 


2.57 


-0.25 


2.00 


-0.63 


2.16 


-0.15 ±0.07 


91 ±6 


£- 


pBoo 


HD 


127665 


K3III 


3.82 ± 


0.12 


21.92 ±0.81 


18.75 


±0.92 


4260 


2.22 


-0.17 


3.58 


0.76 


2.33 


-0.21 ± 0.09 


96 ±9 


re 


i/fUma 


HD 


96833 


K1III 


4.07 ± 


0.03 


22.21 ±0.68 


19.68 


±0.62 


4550 


2.53 


-0.13 


3.01 


0.43 


2.16 


-0.68 ± 0.08 


148 ± 12 


d gian 


y Per 


HD 


9927 


K3III 


3.77 ± 


0.06 


18.76 ±0.74 


21.61 


±0.93 


4380 


2.34 


0.00 


3.57 


0.65 


2.16 


-0.83 ±0.10 


170 ± 16 


24 Per 


HD 


18449 


K2III 


2.10 ± 


0.07 


9.31 ±0.78 


24.21 


±2.19 


4340 


2.37 


-0.19 


4.95 


2.10 


2.33 


-0.73 ±0.18 


155 ±28 




a- Boo 


HD 


124897 


K1.5III 


20.84 ± 0.03 


88.85 ± 0.74 


25.21 


±0.21 


4345 


2.05 


-0.37 


-0.04 


-2.91 


2.16 


-1.01 ±0.05 


201 ± 10 


I 

E 


91Psc 


HD 


8126 


K5III 


2.23 ± 


0.06 


9.49 ± 0.82 


25.30 


±2.29 


4090 


1.93 


-0.17 


5.24 


2.03 


2.51 


-0.58 ±0.19 


135 ±25 


3 


39Cyg 


HD 


194317 


K3III 


3.01 ± 


0.08 


12.77 ± 0.62 


25.35 


± 1.41 


4230 


2.16 


-0.17 


4.44 


1.39 


2.33 


-0.75 ±0.11 


159 ± 18 


re 

s 


11 Lac 


HD 


214868 


K2III 


2.64 ± 


0.05 


10.81 ±0.56 


26.24 


± 1.46 


4440 


2.32 


-0.25 


4.51 


1.67 


2.16 


-1.00 ±0.12 


200 ± 23 


C 


31 Leo 


HD 


87837 


K4III 


3.36 ± 


0.04 


11.89 ±0.72 


30.40 


± 1.88 


4040 


1.81 


-0.02 


4.38 


0.98 


2.51 


-1.14 ±0.14 


227 ±31 




v Boo 


HD 


120477 


K5.5III 


4.75 ± 


0.04 


13.29 ±0.81 


38.44 


±2.37 


3890 


1.55 


-0.23 


4.05 


0.44 


2.51 


-1.44 ±0.14 


299 ± 40 


I 


/?Umi 


HD 


131873 


K4III 


10.09 : 


±0.08 


25.79 ±0.52 


42.06 


±0.91 


4030 


1.83 


-0.29 


2.08 


-1.29 


2.51 


-1.73 ±0.07 


390 ± 25 




aTau 


HD 


29139 


K5III 


20.57 : 


±0.02 


50.09 ± 0.95 


44.13 


±0.84 


3910 


1.59 


-0.34 


0.86 


-2.81 


2.51 


-2.04 ± 0.06 


518 ±32 




y Dra 


HD 


164058 


K5III 


9.90 ± 


0.09 


22.10 ±0.46 


48.15 


± 1.09 


3930 


1.55 


-0.14 


2.23 


-1.16 


2.51 


-1.93 ±0.07 


471 ± 30 




/JCnc 


HD 


69267 


K4III 


5.12 ±0.02 


11.23 ±0.97 


48.96 


±4.23 


4010 


1.71 


-0.24 


3.54 


0.19 


2.51 


-2.05 ±0.19 


526 ± 99 




cfLyn 


HD 


80493 


K7III 


7.45 ± 0.04 


14.69 ±0.81 


54.50 


±3.02 


3880 


1.51 


-0.26 


3.13 


-0.61 


2.51 


-2.32 ±0.13 


673 ± 83 




ySge 


HD 


189319 


M0III 


6.10 ± 


0.02 


11.90 ±0.71 


55.13 ±3.29 


3805 


1.55 


-0.14 


3.53 


-0.16 


2.66 


-2.13 ±0.14 


562 ± 75 




a- Hya 


HD 


81797 


K3II-III 


9.58 ± 


0.08 


18.40 ± 0.78 


55.93 ±2.41 


4120 


1.77 


-0.12 


2.00 


-1.13 


2.33 


-2.48 ±0.10 


780 ± 78 




vUMa 


HD 


98262 


K3III 


4.72 ± 


0.03 


8.88 ±0.64 


57.07 


±4.13 


4070 


1.89 


-0.04 


3.50 


0.28 


2.51 


-2.47 ±0.16 


775 ± 122 





8 LD is the limb darkened angular diameter measured by interferometry, taken from the CHARM2 catalogue compiled by Richichi et al. (2005), except for ySge (Wittkowski et al. 20061, 
S Eri and £ Hya (Thevenin et al. l2005b . e Oph (Mazumdar et al. l2009l and 77 Ser (Merand et al. 120091 . 
The parallax n is taken from the original Hipparcos catalogue (ESA 1997 1. 

r cff , logg and [Fe/H] are taken from Cayrel de Strobel et al. (2001), except for y Sge, 6 Eri, eOph and £ Hya (see text). 
m v and m K were taken from SIMBAD. 

The uncertainty on the K band bolometric correction BC(K) (taken from Houdashelt et al. 2000 except for 7 Sge) is taken uniformly equal to 0.05. 
The absolute bolometric magnitude of the Sun is M bo i(0) = 4.75. 
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4. Convection length scales from the Sun 

The convection length scales can be calibrated using the Sun, 
as is customarily done. We assume L = 3.846 10 33 erg.s~' and 
R Q = 6.9599 10 10 cm and start the solar evolution on the zero 
age main sequence. The calibration in luminosity, radius and 
metal to hydrogen ratio ^saL are achieved to better than 1CT 4 

**-surf 

at the age of 4.6 Gyr for both MLT and CGM convection pre- 
scriptions. We consider the solar abundances and metal repar- 
tition of Asplund et al. (2005) therefore = 1.65 10~ 2 . The 
microscopic diffusion is accounted for following Proffitt and 
Michaud ( 1993 ) It induces a decrease in helium and metal sur- 
face fractions in the course of the solar main sequence. 

In the MLT framework we obtain a m ] to =1.98. The initial 
composition of the calibrated solar model is X Q = 0.7284, 
Y D = 0.2578 and Z D = 1.38 10~ 2 and its solar age sur- 
face composition is X surt = 0.7592 and Y sul f = 0.2282. The 
a m i tQ is compatible with the recent estimate of VandenBerg et 
al. (f2008 ) also using non grey atmospheres as boundary con- 
ditions. Using the MARCS atmosphere models (Gustafsson et 
al. l20031 > the latter authors find a m i t0 =2.Ol . Samadi et al. ((2006 ) 
using the Kurucz atmosphere models (Kurucz 1993, see also 
Heiter et al. 12002 ) find a m i t0 =2.53. This larger value is likely 
due to the fact that Samadi et al. start interpolating between 
the atmosphere and interior gradients at the optical depth t = 4 
whereas we start at t = 1 (see section l53T l. The overadiabaticity 
induced by the atmosphere a — 0.5 extends over a larger region 
in their models and has to be compensated by a more efficient 
convection deeper. We recall that a m \ tG depends on the assumed 
solar composition. For the Grevesse & Sauval ( 1998 ) composi- 
tion (|^ = 2.3 10~ 2 ) calibrated Sun, we find 0^0=1.67. 

In the CGM model framework we obtain a C gmo=0.77. 
The initial composition of the calibrated CGM solar model 
is nearly the same as for the calibrated MLT solar model. 
Z = 1.38 10~ 2 is identical while the hydrogen and helium 
fractions very slightly differ from the MLT case: X D = 0.7280 
and Y = 0.2582. So is the current surface composition with 
X surf = 0.7587 and Y surf = 0.22865. Our Qf C g m © is in close 
agreement with the values of authors who used the same for- 
mulation of the CGM. Bernkopf (1998} finds a cgm0 = 0.82 and 
Samadi et al. (2006) find 0.78. 

Theoretical atmosphere models probably account properly 
for the differential changes in the outer stellar regions with T e ft 
and surface gravity. Yet, in the solar case, empirical atmosphere 
relations seem to reproduce the limb darkening better than the- 
oretical ones (Blackwell et al. 119951 ). Following VandenBerg et 
al. ([2008 ) we combined the advantages of theoretical and em- 
pirical T(t) relations. To that extent the T(t) we used is the in- 
terpolated theoretical ATLAS 12 relations TtheoiyW rescaled by 
the empirical Thm74(t) relation for the solar atmosphere from 
Holweger & Mueller (1974) as suggested by VandenBerg et 
al. (120081 : 

T(T) = T the0 iy(T) + — ^— [T HM 74(t) - T theoryj0 (T)] (1) 
leff,© 

If we introduce such semi empirical atmosphere models the 
solar calibration leads to a slightly larger a cgmO =0.80 for the 
CGM prescription. The initial and final solar surface composi- 



tion change negligibly. We will discuss the influence of such a 
semi empirical atmosphere model in section l5~3l below. 

It is worth giving a word of caution on solar calibrations 
of convective length scales. They are as reliable as are the so- 
lar models. The dynamics of the Sun radiation zone is largely 
uncertain as suggested by helioseismology (Turck-Chieze et 
al. 120041 ). It probably involves angular momentum transport 
through internal waves (Charbonnel & Talon [20031 Mathis et 
al. [2008 ). Moreover the mixing of the radiative interior asso- 
ciated with rotation (as in the tachocline) interacts with diffu- 
sion along the evolution (Brun, Turck-Chieze & Zahn 1999). 
This changes the actual surface composition. In this respect, it 
is interesting to remark that the calibration of the solar mod- 
els using its recent composition determinations have led to he- 
lium content of the convection zone that systematically under- 
estimate the measurement made thanks to seismology (Basu & 
Antia [1995l 

5. Red Giant branch calibrations 

5.1. The cool edge of the RGB 

Depending on its mass, a red giant will not be at the same ef- 
fective temperature for a given luminosity. In the case of open 
or globular clusters all the objects on the RGB approximately 
have the same mass and the RGB is well defined. This is be- 
cause clusters gather stars of the same age. No such assumption 
can be made for a sample that is representative of field stars 
in the solar neighbourhood: stars of different ages (and thus 
masses) are presumably on the RGB we consider. 

It is however possible to set a lower limit on these masses. 
The limit is given by the age of the local Galactic disk and the 
evolutionary time scale of its low mass stars: 

Age: For objects in the slightly subsolar metallicity range 
(-0.25 < [Fe/H] < -0.14) Liu & Chaboyer d2000l ) suggest a 
maximum age of 11.7 + 1.9 Gyr. In this study we mainly con- 
sider models that have reached 1O 3 L on the RGB by that age 
as our RGB stars exhibits [Fe/H] = -0.17 on average. We can 
further set upper limits to the local Galactic disk age: it is prob- 
ably younger than the oldest Galaxy globular clusters whose 
current age estimate is 12.6 Gyr (Krauss & Chaboyer 2003 ) and 
certainly younger than the Universe 13.7 + 0.1 Gyr (Komatsu 
et al. [2009b . 

Initial composition: The average metallicity of the 34 gi- 
ants and 4 subgiants in TableQ]is [Fe/H] ~ -0.17. As is well 
known the first dredge-up occuring near the base of the RGB 
erases most of earlier diffusion effects so that the actual surface 
metal fraction is also the initial one. We set the initial helium 
content at Y = 0.2582 i.e. the amount of our solar calibrated 
models. This choice could be criticized as the stellar nucle- 
osynthesis simultaneously increases the metal and helium frac- 
tions and both of them appear (loosely) correlated (Fernandes 
et al. |1998l >. Furthermore the estimate of the initial solar con- 
tent helium is still a matter of debate as it is affected by the dy- 
namical phenomena in the solar interior (cf final caveat of ® . 
However a discussion on helium would be irrelevant here be- 
cause even a substantial change in its initial fraction has a neg- 
ligible impact on the position of the RGB as we will see shortly. 
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Table 2. The three main different metallicity models. 



J 2.0 



1.5 



[Fe/H] 


Metals 


Hydrogen 




mass fraction 


mass fraction 


0.0 


1.20 10- 2 


0.7297 


-0.17 


8.18 10~ 3 


0.7336 


-0.34 


5.55 10" 3 


0.7362 



i.o r „ i 

3.70 3.65 3.60 3.55 

log(T.„) 

Fig. 2. Distribution of the current stellar sample in the HR di- 
agram. Triangles: lower envelope of the RGB, see text for the 
names of the six corresponding objects. Diamonds: the other 
giants. Star: (5'Tau. 

3.0 '/' ' ' \ 




21og(R/R G ) 

Fig. 3. Distribution of the current stellar sample in the log lu- 
minosity vs log squared radius diagram. The solid line is a lin- 
ear fit to the sample by minimization of^- 2 . Triangles: group 
II, the nine giants having the largest radius with respect to the 
corresponding fitted radius of the sample. Diamonds: the other 
giants. Star: tf'Tau. 



In the bolometric luminosity vs. stellar squared radius diagram 
the lower envelope of the RGB is harder to see: Fig. [3] We 
selected the giants having the largest radii at a given luminos- 
ity in the following manner: first we performed a linear fit to 
the whole set of data. Then we kept only the stars with the ten 
largest differences in radius to the fitted radius at similar lu- 
minosities: cr Gem, 39 Ari, pBoo, 24 Per, 91 Psc, 39 Cyg, 31 
Leo, uBoo, jSUMi and y Dra. As in the case of the classical 
HR diagram 5'Tau was excluded from the sample. ySUMi was 
also excluded for its metallicity is very low ([Fe/H] = -0.29) 
and we refer to the nine remaining stars as 'group IF. Out of 
six stars making the lower envelope of the RGB in the classical 
HR diagram five also belong to the group of the nine most ex- 
panded objects. One could suggest that for simplicity we rely 
on the coolest stars also in the subsequent analysis using the 
radii. We keep groups I and II instead. One purpose of this work 
is to use data on effective temperatures and radii independently: 
the radii are determined through interferometry while the effec- 
tive temperatures are determined through spectroscopy. Within 
the groups, the stars can be further classified depending on the 
metallicity H cr Gem, 6 And, 31 Leo, 39 Ari all are within 0.04 
dex of the solar metallicity. All the other members of groups 
I and II are within 0.06 dex of [Fe/H] = -0.17 (and within 
0.03 dex of this value if v Boo is not considered). Let us now 
compare the data on these objects to the evolutionary tracks of 
models having different compositions, masses and convection 
treatment length scales, first in the case of the MLT and then in 
the case of the CGM model. 



Metals on the opposite have a very strong influence on the ef- 
fective temperatures and radii of RGB stars. Thus in the fol- 
lowing we repeat most of the calculations for the metallicities 
[Fe/H] = -0.17 and [Fe/H] = 0. The fraction Z = 1.20 10~ 2 
of the solar metallicity models very slightly differs from our 
solar calibrated models Z = 1.26 10^ 2 because of the helium 
fraction difference between both sets. In l5.4l we performed ad- 
ditional calculations for [Fe/H] = -0.34 and built models with 
compositions dedicated to the few stars of the sample having 
seismic data. The main compositions investigated are given in 
details in Table[2] 

Between ~30 and ~300 L , six stars clearly define the 
lower envelope of the local Galactic RGB in the classical HR 
diagram: cr Gem, 5 And, p Boo, 91 Psc, 31 Leo and v Boo. 
Hereafter we refer to them as 'group I'. Fig. [2] displays them 
and also shows an object that is far off the general trend. This 
is c^Tau and it will not be considered in the following analysis. 



5.2. The cool edge of the RGB and the MLT 

Table[3]sums up the properties of the different models built us- 
ing the MLT while Fig.|4]and|5]display their RGB evolutionary 
tracks. In the Table, column 1 is the models mass, column 2 the 
metallicity, column 3 the helium mass fraction, column 4 the 
age in Gyr at 10 3 L o , column 5 the a m \ t . Columns 6 and 7 are 
the x 2 goodness of fit between models and data respectively in 
the HR diagram and the radius luminosity diagram. x\ m col- 
umn 6 is based on effective temperatures of stars of group I. 

xt i r rmod_ r robs 

It is defined by X \ = Zjli s[^t^] where T °ff and T ™ 
are respectively the observed effective temperature of an object 
and that of the model corresponding to the same luminosity as 

3 We did not find error bars on metallicities of group I and II stars 
in the catalog of Cayrel de Strobel (2001 1. We just mention that the 
error bars, when present in this catalog, are generally around 0. 1 dex 
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3.70 3.65 3.60 3.55 

log(T.„) 



Fig. 4. RGBs in the HR diagram for models mltl to mlt6 
of Table [3] The solid line is model mltl, the dotted line is 
model mlt2, the dashed line is model mlt3, the dot-dashed line 
is model mlt4, the three-dot dashed line is model mlt5 and 
the long dashed line is model mlt6. Models mlt2 and mlt3 
show mass and helium fraction change with respect to model 
mltl. They are hardly discernable from this model. Model mlt4 
shows the metallicity effect with respect to model mltl. Only 
models mlt5 and mlt6 with lower a m \ t than the previous models 
fit the data satisfactorily. The data are the six giants of 'group 
P defined in the text and Fig. [2] Square symbols are for ob- 
jects within 0.04 dex of solar metallicity while triangles are for 
objects within 0.06 dex of [Fe/H] = -0.17. The error bars in 
log(T e ff) are taken from Mc William (1990), the error bars in 
log luminosity are adapted from Table [TJ 

the object. Cayrel de Strobel et al. (1200 11 1 catalog mostly con- 
siders the T eff from Mc William ([19901 Accordingly our AT* S 
are adapted from this author to AT* S =130 K. x\ in column 7 
is computed from radii of the nine stars of group II according 

to x\ = Efci ^[ Rm A d R 2 R " bs ] 2 - We do not consider the traditional 

obs 

X 2 formula but the reduced^ 2 to allow a comparison of the fit 
when the number of stars is changed. We can see that changing 
the sample of stars from group I to group II hardly changes the 
goodness of the fit. Columns 6 and 7 are based on the the stars 
within 0.06 dex of the metallicity given in column 2 for consis- 
tency between the composition of the models and the observa- 
tions (the nearly solar metallicity objects over groups I and II 
are cr Gem, S And, 31 Leo, 39 Ari). However the^- 2 in paren- 
theses in column 6 and 7 relate the models to all the members 
of group I and II respectively. Finally column 8 is the model 
name and bold fonts indicate modelling inputs changes with 
respect to the (reference) model mltl. 

The following remarks can be done: 

i) The model mltl uses the solar calibrated a m \ t and reaches 
10 3 L G at wl 1.5 Gyr i.e. nearly at the age estimated for the local 
Galactic disk. Its RGB track does not fit the six coolest objects 
within their observational error bars. The models mlt2 and mlt3 
are respectively less massive and helium poorer than the model 
mltl. Model mlt2 is extreme in the sense that it reaches 1O 3 L 
at »13.9 Gyr (or 10 2 L Q at «13.88 Gyr) which is older than 
the current age estimate of the Universe of 13.7 ± 0.12 Gyr 



3.0 




1.0 

1.5 2.0 3.5 3.0 3.5 

31og(R/R Q ) 



Fig. 5. RGBs in a bolometric luminosity vs. stellar squared ra- 
dius diagram for models mltl, mlt4, mlt5 and mlt6 of Table [3] 
The lines conventions are similar as in Fig. [4] Models mlt2 and 
mlt3 are not shown as they could hardly be distinguished from 
the track of model mltl. The data are the nine giants of group 
II defined in text. Like in Fig.|4]square symbols are for objects 
within 0.04 dex of solar metallicity, triangles for objects within 
0.06 dex of [Fe/H] = -0.17. The best fits are obtained for the 
models mlt5 and mlt6. 



(Komatsu et al. 120091 1 . Model mlt3 is also extreme in the sense 
that its helium fraction is the Big Bang Nucleosynthesis one 
(Coc et al. 2004) and evidently cannot be lowered any further. 
Both models mlt2 and mlt3 are in slightly better agreement to 
the data than model mltl and are almost undiscernable from it 
in Fig. [4] They demonstrate that acceptable changes in mass or 
helium fraction can hardly improve the poor fit of model mltl. 

ii) Even though group I exhibits [Fe/H] » -0.1 on aver- 
age, one can distinguish two metallicity subgroups as cr Gem, 
6 And and 3 1 Leo are nearly of solar metallicity while p Boo, 
91 Psc and v Boo are within 0.06 dex of [Fe/H] = -0.17 (see 
Tabl^D- It 1S necessary to investigate the metallicity effects. 
The solar metallicity O.95M model mlt4 provides a better fit 
to the observations than the previous ones: ix&x\ fit to the whole 
group I is 0.95 and it is 0.72 to the subsample of the three so- 
lar metallicity stars. Together with models mlt2 & mlt3, this 
model should also be considered as extreme because of its age 
at 10 3 L o : 13.7 Gyr. The increase in metal fraction lowers the 
luminosity and thus slows down the evolution with respect to 
the same mass model mltl. 

Hi) Models with lower or m i t than the solar value provide 
much better fits to the data. The model mlt5 (a m \ t = 1.68, 
O.95M and [Fe/H] = -0.17) exhibits ^ 2 =0.16. It reaches 
1O 3 L at 11.5 Gyr. The model mlt6 (a m = 1.68, 1.13M Q and 
[Fe/H] = 0) exhibits ;r 2 =0.078. It reaches 1O 3 L at 7.5 Gyr. 
The pattern of the goodness of the fit is similar if one consid- 
ers the luminosity vs. square radius instead of the HR diagram. 
The x\ an d x\ nave the same orders of magnitudes and varia- 
tions for models mltl to mlt4. The two best fits in x\ are a ls° 
obtained for models mlt6 and mlt5 with 0.54 and 1.1 respec- 
tively. In these cases however the ^- 2 s remain larger than the 
x\s. This is mostly due to the small errorbars in square radius 
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M/M 


[Fe/H] 


Y 


Age at 10 3 L o 


<*mlt 


x] 


xl 


Model name 


0.95 


-0.17 


0.2582 


11.5 


1.98 


3.8(3.2) 


4.0(3.8) 


mltl 


0.90 


-0.17 


0.2582 


13.9 


1.98 


3.0(2.5) 


3.2(3.1) 


mlt2 


0.95 


-0.17 


0.2482 


12.3 


1.98 


3.4(2.8) 


3.6(3.4) 


mlt3 


0.95 


0.00 


0.2582 


13.7 


1.98 


0.72 (0.95) 


1.5(1.4) 


mlt4 


0.95 


-0.17 


0.2582 


11.5 


1.68 


0.16(0.11) 


1.1(0.91) 


mlt5 


1.13 


0.00 


0.2582 


7.5 


1.68 


0.078 (0.070) 


0.54(1.3) 


mlt6 



as compared to the errorbars in effective temperatures (See Fig. 
Hand©. 

iv) Recent studies have suggested that the models with a 
solar calibrated a m \ t can properly describe the RGB (Alonso 
et al. 120001 Ferraro et al. [2006t VandenBerg et al. 120(181) . It is 
likely that our different conclusion stems from the systematic 
use of non-grey atmospheres with a low a m \ t = 0.5 parameter. 
Actually, VandenBerg et al. (2008) rely on the MARCS atmo- 
sphere models where they assume the high value ff m i t = 1.80 
of the interior. Alonso et al. (2000) and Ferraro et al. (2006) 
both rely on the Krishna Swamy ( 1966) T - t atmosphere re- 
lation but they do not mention the optical depth Tb at which 
they connect the atmosphere to the interior. In anycase the ab- 
sence of a low efficiency surface convection region is critica- 
ble. Tb ~ 1 would mean that the outermost convective regions 
are handled with the interior ff m i t parameter which is not sup- 
ported by the Balmer lines profiles (see 95.3b - On the other 
hand Tb > 10 would imply that the Krishna-Swamy T - r re- 
lation is extended down in the convective regime, which is un- 
correct as mentioned by Krishna Swamy ( 1966 ). 

5.3. The cool edge of the RGB and the CGM model 

The MLT is widely used in stellar modelling. However the FST 
approach initialy developped by Canuto & Mazzitelli (119911 ) 
is physically more consistent (Mazzitelli 1999). When com- 
pared to observations it seems to provide a better description of 
stellar convection as supported by many recent studies in seis- 
mology (Basu & Antia 1994; Monteiro, Christensen-Dalsgaard 
& Thompson 1996; Samadi et al. 120061 ) and stellar evolution 
(Stothers & Chin [T997l Montalban et al. l200TT ). For this reason 
we find it important to check possible changes of the charac- 
teristic length scale of the FST model. In a more appropriate 
version of this model, the characteristic length scale is the dis- 
tance to the convectively stable region and therefore is not a 
free parameter. This approach is not implemented in our stellar 
evolution code. It requires substantial changes to compute in 
the same iterations both the convection efficiency and distance 
to the convectively stable region. Instead we use the Canuto, 
Goldman & Mazzitelli (fl996l ) version (CGM) of the FST with 
a characteristic length scale crcgmHp as mentioned in § above. 
The precise description of the version of the model including 
the parameters is given in appendix 2. 



Table|4]sums up the properties of the models built using the 
CGM model while Fig. [6] and Q display their RGB evolution- 
ary tracks. Columns conventions are the same in Tables [4] and 
[3] We did not report here on the models exploring the effects 
of a change in mass or helium fraction (corresponding to mod- 
els mlt2 & mlt3 in 15.2b . Similar changes of these parameters 
change the location of the RGB as little as in the case of the 
MLT. 

The following remarks can be done: 

i) Models cgml and mltl are identical except for their at- 
mosphere boundary conditions and convection prescriptions. 
They both suffer the same drawback as the RGB they define 
are too warm with respect to the observations. Model cgml 
RGB is even warmer than that of model mltl. This is a conse- 
quence of the higher efficiency of adiabatic convection in the 
CGM prescription than in the MLT. In the case of deep convec- 
tive envelopes this produces higher effective temperature ob- 
jects (Mazzitelli, D'Antona & Caloi 1995). As can be seen in 
Fig. [6] for larger luminosities (i.e. for larger convective enve- 
lope) the difference in effective temperatures between models 
mltl and cgml becomes larger. 

ii) Having a solar metallicity, the model cgm2 is 0.17 dex 
more metal rich than the model cgml. However in spite of its 
lower effective temperature and larger radius it also fails to fit 
the corresponding observations. As for the mlt models, the x 1 
fits in Table |4] are related to the subgroups of stars exhibiting 
the same metallicity as the models. For the [Fe/H] = models, 
only cr Gem, 5 And and 3 1 Leo are considered within group I 
while only cr Gem, 31 Leo and 39 Ari are considered within 
group II. Models with masses lower than cgm2 and the same 
metallicity could in principle fit the data. However with an age 
over 13.8 Gyr at 1O 3 L , model cgm2 sets a lower limit in mass. 

Hi) Both models cgml and cgm2 have the solar calibrated 
a cgm value. The sub solar a cgm models cgm3 and cgm4 produce 
a better agreement to the selected samples of stars in the two 
considered diagrams (Fig.|6]and|7]( as can be also seen in Table 
g] The model cgm3 has [Fe/H] = -0.17 and 0.95 M while 
cgm4 has [Fe/H] = and 1.17 M . Thus models cgm3 and 
cgm4 correspond respectively to models mlt5 and mlt6 in terms 
of masses and metallicities. 

iv) In an attempt to investigate the atmosphere modelling 
effects, the model cgm5 was built with the semi empirical ap- 
proach for the atmosphere calculation described in section [4] 
Accordingly its a cgm is 0.80 which is the solar convective scale 
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Table 4. Characteristics of the RGB models using the CGM. 



M/M Q 


[Fe/H] 


Y 


Age at 10 3 L o 


^cgm 


x\ 


x\ 


Model name 


0.95 


-0.17 


0.2582 


11.6 


0.77 


6.1(5.1) 


7.1(6.2) 


cgml 


0.95 


0.00 


0.2582 


13.8 


0.77 


1.1(1.6) 


1.8(2.3) 


cgm2 


0.95 


-0.17 


0.2582 


11.8 


0.62 


0.21(0.11) 


0.70(0.57) 


cgm3 


1.17 


0.00 


0.2582 


6.9 


0.62 


0.057 (0.077) 


0.40(0.76) 


cgm4 


0.95t 


-0.17 


0.2582 


11.5 


0.80 


7.4(6.1) 


9.1(7.8) 


cgm5 



Bold fonts indicate inputs changes with respect to the (reference) model cgml. The model cgm5 mentioned by f is built with the semi 
empirical atmosphere described in section|4] 



when the semi empirical atmosphere are used. The cgm5 RGB 
lies very close to the cgml RGB and is therefore too warm to fit 
the lower envelope of the observations. Models cgml and cgm5 
illustrate that a change to semi empirical atmospheres does not 
reduce the gap between calculations and observations. 

5.4. The global RGB and the seismic targets 

We now intend to use the rest of the data we have about the 
sample of 38 stars. First, we check how the stars are distributed 
between the evolutionary tracks of different masses. We con- 
sider the two surface convection prescriptions with convective 
length scales calibrated on the Sun or on the lower envelope of 
the current RGB in the preceding section. 

Assuming that the stellar mass distribution of our sample is 
the present day Galactic field stellar mass function (hereafter 
PDMF), the question we can investigate is whether the distri- 
bution of the sample between the evolution tracks of different 
masses agrees with the PDMF. The use of the PDMF here re- 
lies on two assumptions. First, it assumes that the mass loss 
of the preceding main sequence and giant branch evolution is 
negligible. There is no indication of significant mass loss be- 
fore the RGB tip. When using the recent formula advocated in 
Catelan (2000) we find that the 1.5M model looses less than 
0.005M o before reaching 8 1O 2 L . For lar ger masses, the mass 
loss is lower. Second, it assumes that the stars of our sample are 
a well mixed Galactic field population that formed in a variety 
of star forming regions and environments. Should a substantial 
fraction of the stars in our sample be formed in a single proto- 
cluster star forming region, deviations, both in the slope and in 
the characteristic mass, between the PDMF and that of our sam- 
ple could be expected (Dib et al. 2010). The observed PDMF of 
stars in the local Galactic field is well fitted by a multi-exponent 
power law functional form (Kroupa 2002 and 2007]!: 

— = k(0.5/0.08r L3 (M/0.5r 2 ' 3 if 0.5 < M/M < 1.0 (2) 
dM 

— = k(O.5/O.O8r L3 (l.O/O.5r 2 - 3 (Mr 2 - 7 ifM/M > 1. (3) 
dM 

In Eq. |2] and [3] k=0.877 and N is the density of stars of 
mass M (in solar mass units). The corresponding repartition 
below between and above 1.5 and 2.5 M is given in the last 



column of Table[5] We name it 'truncated' PDMF as we derive 
the distribution using the previous power laws and making the 
hypothesis that no star has a mass lower than O.95M which 
we found to be the mass of the objects making the lower enve- 
lope of the RGB ( 35.21 i. If we had considered a lower trunca- 
ture mass such as O.9OM (corresponding to an unrealistic age 
for the oldest stars of the sample) the distribution would hardly 
have been affected. As is well known, the position of the RGB 
strongly depends on the metallicity. For instance, for the 1 .5M 
CGM models at 1O 2 L , a change from [Fe/H] = -0.34 to 
[Fe/H] = lowers the T e ff by 240 K (and respectively increases 
the radius from 15.8 to 17.7 R ). Thus the distribution of stars 
in regards to evolutionary tracks is made with models nearly 
having the same metallicity as the observed stars. We compare 
data and models for three metallicities. The [Fe/H] = -0.34 
tracks are used for objects exhibiting [Fe/H] < -0.25. The 
[Fe/H] = -0.17 tracks are used to estimate the distribution of 
objects having -0.25 < [Fe/H] < -0.08. Finaly the [Fe/H] = 
tracks are used when the observation gives [Fe/H] > -0.08. 

The repartition between the tracks has been drawn from the 
luminosity vs. stellar squared radius diagram (see Fig. [8}. It is 
reported in Table [5] for the different assumptions on the con- 
vective length scales. Column 2 shows the distribution for stars 
with a m \ t = 1.98, column 3 for stars with a m \ t = 1.68, col- 
umn 4 for stars with a cgm = 0.77 and column 5 for stars with 
ttcgm = 0.62. For a given mass, composition and luminosity, the 
larger the characteristic convective length scale the warmer the 
model. Being warmer, the solar calibrated models {a m \ t = 1.98 
or a cgm = 0.77) suggest a local RGB strongly biased toward 
lower mass stars. The models calibrated on the lower enve- 
lope of the RGB {a m \ t = 1.68 or a cgm = 0.62) are in better 
agreement with the PDMF. This confirms the slightly less effi- 
cient CGM convection in the RGB regime found in section [531 
However this trend is merely indicative. The small sample of 
stars in each mass bins induces large statistical uncertainties. 
Moreover the RGB tracks of stars with convective cores on the 
main sequence depend on the amount of overshooting during 
this phase. This amount of core overshooting is not a well es- 
tablished quantity. The larger the overshoot the more massive 
the helium cores and the smaller the radius at a given lumi- 
nosity. According to the mass dependence of core overshoot- 
ing of Claret 2007 (see Fig. 13 of this author), we have taken 
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Mass range 


ffmlt 


Omit 


Q'cgm 


^cgm 


Truncated 




1.98 


1.68 


0.77 


0.62 


PDMF 


M < 1.5M 


76 ± 4% 


54 ± 13% 


84 ± 15% 


49 ± 11% 


54% 


1.5M < M < 2.5M 


24 ± 8% 


41 ± 11% 


16 ± 7% 


40 ± 10% 


27% 


M > 2.5M S 


0% 


5 ±4% 


0% 


11+5% 


19% 



a m = O.lHp and a ov = 0.2H p in 1.5M and 2.5M Q models 
respectively. 

We dedicated a special attention to 6 En, £Hya and eOph 
as asteroseismology has allowed an accurate determination 
of their masses. Our 5Eri models have a mass of 1.215 M 
(Thevenin et al. 2005) and X=0.7256 to account for the star's 
metallicity [Fe/H] = 0.13. The eOph models have a mass of 
1.85 M Q (Mazumdar et al. [20091 and X=0.7326 ([Fe/H] = 
-0.12). Finaly the £Hya models have a mass of 2.65 M 
(Thevenin et al. 2005 ) and a hydrogen mass fraction X=0.7308 
([Fe/H] = -0.04). All those models have the same helium 
mass fraction Y=0.2582. We assume a convective core over- 
shooting of 0.1 Hp in <SEri and 0.2 H p in £Hya and eOph. 
Fig(9] shows that the position of the 5 En models weakly de- 
pends on the chosen convection length scale. The a cgm = 0.77 
and a m it = 1 .98 (solar calibrated) models cannot be ruled out 
as they fit the data in the errorbars. However the a cgm = 0.62 
and a m ii = 1 .68 models converge right to the observed radius 
at the observed luminosity! For £Hya only the a cgm = 0.62 
and Ofmit = 1.68 models fit the observations within the error- 
bars: Fig[10] On the contrary, the eOph tracks with the low 
Ccgm = 0.62 does not fit them. It is a very interesting point. 
The explanation might be that e Oph is not on its first ascent of 
the RGB but a 'red clump' star i. e. in the helium core burning 
phase. There are several clues of this. First eOph belongs to 
the tight group of six stars that can be seen at log(L/L ) « 1 .7 
on Figj3]and that strongly suggest the red clump. The second 
clue comes from Mazumdar et al. ( 120091 ) models. They spend 
20 times more time in e Oph errobox in the HR diagram while 
on the helium core burning stage than while on the the first as- 
cent of the RGB which means it is much more likely that we 
observe an helium burning e OphQ. Finaly, except for its mass, 
e Oph seems identical to £ Hya which is likely on its first ascent 
of the RGB. Their metallicities only differ by 0.08 dex. Their 
luminosities (5% difference), radii (2% difference) and effec- 
tive temperatures (80K difference) are the same within error- 
bars but the masses strongly differ. Such a difference induces a 
difference in models radii of « 15% at « 6OL if both stars are 
on the RGB and if the same convective length scale is taken. 
Contrary to e Oph, £ Hya is probably not in the helium burning 
stage. We found the minimum luminosity of the 2.65 M mod- 
els during the core helium burning (eAGB) around 80 L . This 
value, in agreement with stellar grid calculations (Schaller et 

4 Our stellar evolution code cannot pass the so-called helium flash 
and follow the subsequent evolution of stars lighter than w 2.1M G like 
6 Oph during their helium core burning phase. 



al l 1992| l, stands well outside the luminosity errorbar for £ Hya. 
Finaly, leaving apart the case of e Oph, the seismic data com- 
bined to radii measurements suggest a drop in the convective 
length scale confirming the previous results. The seismic tar- 
gets are in small number at the moment and Fig|9land[T0lmostlv 
illustrate that combining seismologically determined masses to 
interferometric radii offers very sensitive tests of the surface 
convection efficiency. In the near future seismology should al- 
low many more accurate mass and radius determinations (Basu, 
Chaplin & Else worm [2010]). 

6. Discussion 

This work aims at constraining the outer convection prescrip- 
tions in low mass red giants and subgiants by using their abso- 
lute luminosities, effective temperatures, interferometric mea- 
surements of their radii and seismic data. The observational 
sample is made out of 38 Galactic disk nearby stars. It was 
selected on the basis of interferometric radii measured with 
a better than 10 percent accuracy. The average metallicity of 
the sample is [Fe/H] = -0.17. There are small but signifi- 
cant metallicity differences between the stars. Age and mass 
differences are expected as well. From the modelling point of 
view, we used a modified version of the secular stellar evo- 
lution code CESAM. A special care was devoted to the at- 
mosphere boundary conditions. We computed two grids of 
non grey atmospheres surrounding the expected surface grav- 
ities and effective temperatures. The first grid, based on the 
PHOENIX/ID calculations, relies on the mixing length theory 
for the treatment of convection while the second grid, based 
on the ATLAS 12 calculations, relies on a modified version of 
the full spectrum of turbulence model (see appendix 2). We 
considered the boundary conditions to the interior model at the 
Rosseland optical depth 20. In the regime of surface conditions 
encountered in red giants, superadiabatic convection extends 
above and below this limit. The same model of convection was 
used for the atmosphere and the interior but different convec- 
tion length scales were adopted in the two regions. A procedure 
of linear interpolation of the thermal gradient with the optical 
depth allows a smooth transition of it between them. Other im- 
portant changes to CESAM were made to enable modeling of 
late stages of stellar evolution (see appendix 1). We chose the 
Asplund et al. (2005) solar repartition for metals in the opacity 
tables. The total metallicity was varied to account for solar or 
slightly sub solar metal content in the sample and the variations 
between the stars in dedicated cases. 
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3.70 3.65 3.60 3.55 

log(T.„) 

Fig. 6. RGBs in the HR diagram for models cgml to cgm5 of 
Table [4] The track of model mltl of Table [3] is also given for 
comparison. The solid line is model cgml, the dashed line is 
model mltl, the dot-dashed line is model cgm2, the three-dot- 
dashed line is model cgm3 and the long dashed line is model 
cgm4. The dotted line is model cgm5. The data are the six gi- 
ants of group I defined in the text and Fig. [2] Following Fig. [4] 
the squares and triangles respectively stand for solar and subso- 
lar metallicity objects. Only models cgm3 and cgm4 with lower 
a cgm than the other models fit the data satisfactorily. 
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1.5 2.0 2.5 3.0 3.5 

21og(R/R s ) 

Fig. 7. RGBs in a bolometric luminosity vs. stellar squared ra- 
dius for models cgml to cgm4 of Table|4] The lines conventions 
are similar as in Fig. [6] The data are the nine giants of group II 
defined in the text (see also Fig. [3). The data symbols are the 
same as in Fig. [5] Triangle symbols are for objects within 0.04 
dex of solar metallicity, triangles for objects within 0.06 dex of 
[Fe/H] = -0.17. 

We proceeded in three steps. First we built calibrated so- 
lar models for the two prescriptions of convection: the MLT 
and the CGM version of the full spectrum of turbulence. Then 
we built RGB models for different masses, metallicities and 
initial helium fractions. The comparison between models and 
data was made in the classical HR diagram and in the squared 
radius vs. luminosity diagram. In the third step we investigated 
the mass distribution on the red giant branch. To that extent 
we made the assumption that the mass distribution of the data 
set follows the present day Galactic field stellar mass function. 




0.5 - ♦ • 



o.o h 1 

12 3 4 

21og(R/R Q ) 

Fig. 8. Exemple of RGBs in the bolometric luminosity vs. stel- 
lar squared radius diagram for models using the CGM convec- 
tion model and having [Fe/H] = 0: 1.5 M solid line and 2.5 
M dashed line. The evolutionary paths are shown from the 
zero age main sequence. For the sake of clarity only the stars 
with [Fe/H] > -0.08 are displayed. 
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0.60 0.65 0.70 0.75 0.80 
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Fig. 9. 5 Eri evolution tracks for a cgm = 0.62 (solid line) 
ffcgm = 0.77 (dotted line), a m \ t - 1.68 (dot-dashed line) and 
ffmit = 1-98 (dotted line). All models match the observed ra- 
dius at the observed luminosity within the errorbars. However 
the lower convective length scales models a cgm = 0.62 and 
£*mit - 1 -68 better fit the data. 

Finaly we use of the asteroseismically determined masses in 
the few cases where they are available. 

Our conclusions are the following: 

i) Lower envelope of the RGB 

When the mixing length theory is used with the solar con- 
vective length scale (A = a m i t H p and a m \ t = 1.98), the model 
with 0.95M Q & [Fe/H] = -0.17 ascends the RGB at * 11.5 
Gyr i.e. the estimated age of the local Galaxy disk. This model 
poorly fits the lower envelope of the RGB in effective tem- 
peratures or radii. The situation is not improved by lowering 
of the mass or the helium fraction unless we consider models 
older than the Universe or having an helium fraction below the 
Big Bang Nucleosynthesis value. Considering a solar metal- 
licity model of 0.95M o improves the fit but at the expense of 
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1.8 1.9 2.0 2.1 2.2 

21og(R/R G ) 

Fig. 10. Solid lines show the evolution tracks for a cgm = 0.62 
and a m \ t = 1.68 and a mass of 2.65M . The dotted lines cor- 
respond to the same mass but a cgm = 0.77 and a m \ t = 1.98. 
The solar calibrated models for convection do not fit the obser- 
vations. Dashed line and dot-dashed lines are 1.85M models 
respectively for a cgm = 0.77 and a cgm = 0.62. The lower radius 
object is £ Hya and the larger e Oph. 



an age (13.7 Gyr) hardly compatible with current constraints 
on the age of the Galaxy. Models with lower convective length 
scale {a m \ t = 1.68) much better fit the observations both in ef- 
fective temperatures and radii. Those models suggest two age 
and mass components for the lower envelope of the local RGB. 
The 0.95M o & [Fe/H] = -0.17 model reaches 10 3 L G at 11.5 
Gyr and the 1.1 3M & [Fe/H] = model reaches this lumi- 
nosity at 7.5 Gyr. Interestingly this result is in agreement with 
the work of Liu & Chaboyer ( 120001 ) who estimate that the age 
of field stars showing [Fe/H] = -0.18 is 11.7 + 1.9 Gyr and 
the age of field stars with [Fe/H] = is 7.5 +1.7 Gyr (see the 
Table 5 of these authors). 

When the full spectrum of turbulence model is used with 
the solar convective length scale (A = ffcgmHp and a cgm = 
0.77), the O.95M & [Fe/H] = -0.17 model is also too warm 
with respect to the lower envelope of the red giant branch. 
Correspondingly it gives too small radii when compared to 
the observed ones. This feature is expected as when the con- 
vection is efficient (e.g. in giant stars envelopes), the CGM 
model convective flux is roughly ten times larger than with the 
MLT one (Canuto, Goldman & Mazzitelli I1996K Eventually 
this produces higher effective temperatures. One needs a con- 
vectively less efficient model a cgm - 0.62, with O.95M & 
[Fe/H] = -0.17 to provide a good fit whether the surface tem- 
perature or the radius is used. Like in the MLT case the CGM 
models suggest two ages for the cool edge of the local RGB, 
the low metallicity component ([Fe/H] = -0.17) being 11.8 
Gyr and the high metallicity one ([Fe/H] = 0) being 6.9 Gyr. 
We did not test the more appropriate version of the full spec- 
trum of turbulence with the convective length scale A equal to 
the distance z to the convective boundary. 

ii) Global RGB and seismic targets 

We studied the mass repartition on the local RGB using 1.5 
and 2.5 M models. We found that models having the solar 



calibrated convective lengths predict a mass distribution on the 
RGB that seems incompatible with the local stellar mass dis- 
tribution function. The lower a m \ t — 1 .68 and a cgm = 0.62 are 
in better but not complete agreement with the expected mass 
distribution. Stars with seismic mass determination offer a final 
way of probing surface convection efficiency. Only two stars 
of the sample could be used that way. The models dedicated to 
them are in excellent agreement with the constraints on lumi- 
nosities and radii. Confirming the previous results, they support 
a slightly lower than solar surface convection efficiency on the 
RGB. 

To describe the RGB the convective length scale needs to 
be slightly decreased with respect to the Sun in the versions of 
the MLT and the CGM we used. This result is suggested in- 
dependently by effective temperatures, radii measurements and 
the few constraints from asteroseismology. The drop is in dis- 
agreement with previous analyses on the red giant branch. We 
interpret this point as a consequence of our systematic use of 
atmosphere models with low efficiency convection. In surface 
layers, the specific entropy increases inwards from the point 
where the convection sets on until the region where the convec- 
tion becomes adiabatic. In a solar model 35% of this increase 
occurs above the optical depth r = 20. But on the RGB typi- 
cally 15% of the increase occurs above t — 20 (estimate from 
model cgm4 at 100 L ). Going from solar to RGB surface con- 
ditions means that the superadiabatic convection -and therefore 
the radius of the star- depends less and less on the atmosphere. 
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Appendix 1 : new CESAM integration variables 

It is worth mentioning an important change in the integration 
scheme of the stellar evolution code we use. The hydrostatic 
equilibrium equation (^ = -^?r) and the radius-mass equa- 
tion (jjgj = 4^2) have singularities at the center of a star. 
For this reason and also because the CESAM code uses piece- 
wise polynomials to describe the stellar structure (Morel 1997 ), 
the integration variables chosen initialy are m 2 ^ 3 , £ = r 2 and 
A = £ 2 I 3 where m, r and I are respectively the local lagrangian 
mass coordinate, radius and luminosity. This choice avoids the 
singularities and makes the calculations stable. However it pre- 
vents the calculation of the stellar structure in case the lumi- 
nosity becomes negative. 

Such a situation arises in the helium cores of low mass stars 
when they ascend the red giant branch. Then the core is devoid 
of any energy source (but its own slow contraction) but effi- 
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ciently looses energy because of the neutrinos. At some point 
the hottest region in the star becomes a spherical shell which 
generates an inward radiative energy flux below i.e. a nega- 
tive luminosity. We changed the integration variables of the 
CESAM code to m 1 ^ 4 , £ = r 2 and A = t. This choice leaves no 
singularities is stable and allows the calculation of the structure 
even though the luminosity becomes negative. To our knowl- 
edge this choice of integration variables has not been made with 
the CESAM code in previous calculations. 



Appendix 2: the current FST convection treatment 

There are several versions of the MET as well as of the FST 
model. Whenever such treatments of convection are used it is 
important to describe them and to provide their parameters in 
order to make precise comparisons to other works possible. 
The current MET convection treatment is described in detail 
in Piau et al. (2005) whereas the current FST version relies 
on the Canuto, Goldman & Mazzitelli (1 19961 1 equations. The 
same equations are used both for the atmosphere and the inte- 
rior convection. The atmosphere models were computed with 
the ATLAS 12 code (Castelli 120051) and especially for the ac- 
tual solar surface composition i.e. X=0.7392, Z=0.0122 and 
the metal repartition advocated by Asplund et al. (2005 ). The 
convective flux F conv is given according to: 



7 conv - K rac ]TH p (V ■ 



where K la d = is the radiative conductivity and S is the 
convective efficiency: S = Ra x Pr oc A 4 . Ra and Pr are re- 
spectively the Rayleigh and Prandtl numbers of the convective 
flow while A is its characteristic length scale. All other symbols 
keep their traditional meanings. 

The <1>(S) function is the ratio of convective to radiative 
conductivity and is given by: 



Vad)O(S) 



where 



and 



<D(S) = F,(S)F 2 (S) 



Fi(S) = (K /1.5) 3 aS k ((l + bS) m - l) n 



F 2 (S) = 1 + 



eS r 



cS p 



1 + dSi 1 + fS< 



with K Q =1.7, the Kolmogorov constant, a=10.8654, 
b=0.00489073, k=0.149888, m=0.189238, n=1.85011, 
c=0.0108071, p=0.72, d=0.00301208, q=0.92, 
e=0.000334441,r=1.2,f=0.000125,t=1.5. 

The equations used here are similar to those of Heiter et 
al. ((2002). However unlike those authors we do not adopt A = z 
where z is the distance to the boundary between radiatively sta- 
ble and unstable regions. It would require substantial changes 
of our code to solve the stellar structure equations with A = z. 
The reason is that this distance z is not known before the 
equations have been solved and makes the problem non local. 
Instead A = a C gmH p is a local quantity which leads to a simpler 
integration scheme. 
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